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Abstract 

We consider the axial compression of a thin sheet wrapped around a rigid cylindrical substrate. 
In contrast to the wrinkling-to-fold transitions exhibited in similar systems, we find that the sheet 
always buckles into a single symmetric fold, while periodic solutions are unstable. Upon further 
compression, the solution breaks symmetry and stabilizes into a recumbent fold. Using linear 
analysis and numerics, we theoretically predict the buckling force and energy as a function of the 
compressive displacement. We compare our theory to experiments employing cylindrical neoprene 
sheets and find remarkably good agreement. 

PACS numbers: 46.32.+X, 46.70.-p, 68.60.Bs 
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I. INTRODUCTION 


When you roll up your sleeves to get some work done, you will not pay attention to the 
intricate folding patterns which form around your arms. However, these patterns are not 
only interesting for graphics designers but of eminent importance for biology and technology 
alike iDia- one can hnd them in the twinkling of an eye |3] as well as in the development 
of organs such as the brain 015], the intestine [HIE], or the kidney [8]. Technological 
applications include structures for optics [9] or microfluidics [10] to name just a few. 

One common theme in these examples is that they consist of coupled layers which un¬ 
dergo morphological changes in response to external or internal constraints such as a simple 
compression or volumetric growth. The materials involved range from swelling hydrogels 
mm to supported graphene mm and many others. A well-studied setup in this con¬ 
text consists of a stiff membrane attached to a flat elastic or fluid bulk material HMD. 
The interplay between the bending of the sheet and the response of the bulk leads to the 
formation of wrinkles when the sheet is compressed uniaxially. Beyond a critical compres¬ 
sion the wrinkles vanish and localized folds appear in the sheet. Interestingly, the shape of 
the sheet on the fluid can be found analytically [T8lj2n] as long as the sheet does not touch 
itself [2T] . 

When the substrate is not flat, translational invariance is broken and a whole plethora 
of folding patterns can be found In this article we study a particular type 

of such a system in a cylindrical geometry. An elastic cylindrical membrane is wrapped 
around a solid cylinder of same radius and compressed parallel to the axis of symmetry. 
This simple system is of potential relevance for situations as diverse as intestinal inversion 
m, the folding of your sleeve (see Fig. 1), or even the neck of hidden-necked turtles [3]. In 
contrast to the aforementioned flat system we allow the membrane to stretch azimuthally 
to accommodate to the external stress. Without the solid cylinder constraint the membrane 
would behave like a hollow cylindrical shell whose mechanics has been studied extensively 
in the literature m ESj: when compressed the shell develops regular patterns, such as 
periodic, axisymmetric undulations 1271 or the trapezoidal patterns found by Yoshimura in 
the 1950s [28]. As is easily conhrmed by compressing an empty can of soda, axisymmetric 
modes of deformation are typically unstable. As we will see below, the behavior becomes 
fundamentally different when the shell enwraps a solid cylinder of same size. 
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Similar to a ruck in a rug [29H32] we will consider the case in which the sheet can slide 
on the cylinder. The only coupling between the sheet and the substrate is via the hard 
cylinder constraint. There is no elastic response between the two as is typically the case 
for cylindrical core-shell materials [33]. To simplify the theoretical treatment we assume 
that the sheet is unstretchable in the axial direction. This assumption will be validated by 
experiments with neoprene sheets and hnite element simulations. 

We start with the presentation of our model in Sec. [IT] The resulting shape equations are 
axisymmetric and can be linearized and solved in lowest order of the axial compression as 


shown in Sec. Ill In Sec. IV solutions of the full nonlinear system are found numerically 


and compared to experiments with neoprene sheets. 


II. MODEL 


We consider a cylindrical elastic sheet of thickness h which enwraps a cylinder of radius 


Rq (see Fig. 1(b)). The axis of symmetry is oriented along the X axis whereas we use R 
as the variable for the radial displacement. When the sheet is compressed with a hxed 
displacement AX, it buckles out of its reference conhguration and forms a fold which tips 


over for large compressions (see Fig. 1(c)). In the following we will use the angle-arc length 
parametrization which describes the shape of the axisymmetric sheet with the help of the 
tangent angle 'll) as a. function of arc length S. 

The total elastic energy is given by the sum of stretching and bending contributions. We 
use the linear-strain model described in the appendix with the elastic energy 


E = 


Yhn 
1 — z/2 
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+ 12 
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Rq 
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1 
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1 
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( 1 ) 


We are interested in the case where Ro/h S> 1. The stretching term YR ^o) . clearly 
diverges as H/Rq approaches 0, and similarly does the hrst bending contribution from xp''^. 
The second bending term oc h^Ro{^^^ — ;^)^, however, vanishes for large h/Ro since costp 
is bounded. The last term, which looks like a Gaussian curvature, but is not due to the 
integration over the reference conhguration, can be integrated to a term proportional to 
— sin 'ip E'lp and is constant due to the hxed angle '0 = 0 at the boundaries for all conhgura- 
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(a) 



(b) 


(c) 


FIG. 1: Top: (a) Rolling up the sleeve. Bottom: The studied system (b) before and (c) after the 
tipping point (overlay of a photograph of the neoprene sheet used in the experiments with a sketch 
of the variables of the theory). 


tions. Therefore, we do not need to take it into account either. This is true even for small 
Ro/h: Using linear strains, the problem is independent of the material’s Poisson ratio. We 
are left with the simplihed model 
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We dehne k = and introduce some additional variable rescaling 
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in order to write the energy functional as 


E 


e : = 


2'kRqk Jo 


ds C 


ds 


-(-0'^ + 12p^) + Xp {p' - sinj;) + {x' - cosj;) 


(4) 


where the dash denotes derivatives with respect to s. The Lagrange multiplier functions 
and Xp couple the Cartesian coordinates x and p to ip. The conjugate momenta are 


- — - I' - \ \ 

Pp r ) Px Xx , Pp Xp . 

We switch to a Hamiltonian formulation with the Hamiltonian 


(5) 


= tp'pp + x'px + p'pp - C 
= ^ ~ + Pp sin -0 + Px cos -0 , 


from which we obtain the Hamilton equations: 


X = 


Pp 


Px = 


m 

dip 

on 


" dx 

dn 


Po = 
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( 6 ) 
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dpp ’ 

(7a) 

dH 

= cosip , 

^Px 

(7b) 

dn . , 

— = sm V’ , 

OPp 

(7c) 

= Px sin Ip - Pp cosip , 

(7d) 

= 0, 

(7e) 

= 12p . 

(7f) 


This set of equations strongly resembles the geometrically nonlinear Euler beam model, 
which reads |3l] 

tp” = F sin p) — Fn cos ip . 


Taking the derivative of Eq. (7a) and combining it with Eq. (7d), we get a ‘modihed’ Euler 
beam with 

tp" = Px sin Ip — Pp cos tp 
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in scaled variables. The essential difference is that the resulting ’normal’ force component 
Pp is now a function of s. The Lagrange multiplier p^ is a constant along the contour due 


to Eq. (7e) and is directly related to the external horizontal force F : 

F 


Px = 


2ttRk 


y/hR=: -f , 


( 8 ) 


where / is positive when the sheet is compressed. 

To find the shape of a single fold, the Hamilton equations ([^ have to be solved with the 
appropriate boundary conditions: 


p(0) = p{sc) = 0 , (9a) 

-0(0) = ^/J{sc) = 0 , (9b) 

'0'(O) = ^|J'{sc) = 0 , (9c) 

n = -f. (9d) 


Eq. (9b) takes into account that the membrane must not have kinks: at the contact point 
s = 0 the membrane detaches from the cylinder and -0 equals 0; at s = Sc, we have p{sc) = 0 


and the profile is horizontal again. The contact curvature condition (9c) results from an 
energy balance at the contact point |35l|36], whereas the Hamiltonian (9d) is a constant due 
to the fact that we do not fix the total arc length Sc- 


III. LINEARIZATION 

Owing to the strong nonlinearities, we are unaware of any analytical, closed-form solution 
of Eqns. (j^. Nonetheless, certain unknowns such as the buckling force and the fold length 
can already be obtained in good approximation using linearized equations. To derive them, 
we start from Eqns. ([^ of a single fold and eliminate p by taking the derivative of Eq. ([7f[). 
This yields 


Pp = 12 sin -0 , 

ijj'' = —f sin 0 — Pp cos -0 . 


( 10 ) 

( 11 ) 


We define 



( 12 ) 
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and reparametrize the equations above using s := 2s/sc — 1 to obtain 


Pp = 3Sc sin?/^ 


r 


1 
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cos V’ ds 


'-1 


—/sin-i/) 
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Pp cos -0 , 


( 13 ) 

( 14 ) 

( 15 ) 


where / = fs^/A is positive since the shell is axially compressed. Dashes denote derivatives 
with respect to s from now on. We are interested in solutions where e <C 1. Therefore, we 
use the following expansion for pp, ^jJ, F and sp. 


^ + 0{e^) , 

( 16 ) 

Pp=P^P +Pf + 0{€^) , 

( 17 ) 


( 18 ) 

Sc = 5^°) + + 0{e^) . 

( 19 ) 


Plugging the expansions into Eqns. (13)/(14) and collecting terms of linear order in e 
gives 


p<‘'" = 3(4»>)>"> , 

^(1)" = _/(0)^(l) _ . 

Combining both equations into one 4th order equation, we hnally get 

o/ (0)\4 


the general solution of which is given by 


= cieVM + c2e-V"+^>C3eV"-^' + 


046 


with A| = 


-/»±C(/(o))2-3(s®)4 


Note that is imaginary since > 0. One directly obtains: 


/>"> = |A5_| + |A= I and s™ = 2{J\\l\\\l\ll2 . 


( 20 ) 

( 21 ) 

( 22 ) 

( 23 ) 

( 24 ) 

( 25 ) 
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For small to moderate displacements, we expect symmetric folds. We thus assume 
"0 to be asymmetric in s —>■ —s and construct asymmetric solutions = 

^ (V'gen(s) - V’gen(-'S) ), leading to 


(26) 


= Oi sin + 02 sin , 


where the coefficients ai = i (ci — C2) and 02 = i (03 — C4) are real. By restricting ourselves to 
these solutions, it suffices to satisfy boundary conditions at one side. Enforcing = 0, 

we can determine Oi as 

sinfv^ 


Ol — —0,2' 


sm 




(27) 


leading to 




-CSC 


yhh) sin (\/|hi) sin (^/lyjs) + sin (^/lEjs 


(28) 


We hrst consider solutions where the curvature = 0 at s = 1. This leads to the condition 


yiTi cot - yTh not = 0 


(29) 


We recall that the Hamiltonian of the folded sheet is conserved along s. This provides an 
additional equation to determine admissible values of -^/lA^I. In rescaled variables Eq. ( 6 ) 
reads 


H = 




n 


c 2 , -I / 

bp +PpSmip - -cosip, 


(30) 


where p is the solution to p' = ^ sin-^. From the boundary condition (9c) we hnd "H = — ^ 


for — 1 < s < 1. In particular in the middle of the fold (s = 0) we thus obtain for each order 
in e 


f sinV-ds ) = 0 . 


'-1 


Fig. H shows the roots of Eq. ( [^ (solid black) and ( [3l| ) (dashed black) as a function 
of and the buckling force Crossing points correspond to solutions satisfying both 
conditions. The dotted red lines in Fig. [^denote singularities of Eq. (29), which excludes, 
e.g., the point (si°\/*^°^) ~ (4.8,8) from the set of solutions. Upon inspection, we observe 
that solutions are found for -^/lA+l = ^ and a/|A^| = ^ where m > n > 0 are impair 


integers. According to Eq. (25), the choice m = 3 and n = 1 yields the lowest possible 






























FIG. 2: Roots oi % + (dashed) and (solid black) as a function of and 

Crossings between both roots correspond to nontrivial solutions of the linearized fold equation. 
Dotted red lines denote singularities of The solution with the lowest force is found for 

(sc^\ ~ (2.92,11.55). The inset shows the buckling force as a function of the rescaled curvature 

boundary condition (1) = ci/a 2 . The force reaches a minimum for ci = 0. The dotted line 
shows the buckling force for unphysical positive values of ci/a 2 , where the fold would buckle 
inwards. 


buckling force and thus corresponds to the physical solution with 


m ^ ^ 
^ 2 


24.67, 


= 1 V 


1/4 


TT 2.92 


or = 20/v^ ~ 11.55. The angle at linear order follows as 

Svr, 
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= 02 1 sin —s + sin —s 
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(32) 


(33) 


It should be noted that the amplitude 02 is still undetermined at this point. It can be 
found by prescribing the displacement e^: Integrating the displacement constraint equation 


Eq. (15) and retaining terms up to order yields 


1 - 6^ ^ - r 

1 

1 
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bO 

_1 

2 7-1 
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1 ~ 1 ^2 

ds = i-Y. 


(34) 


and thus 02 = — \/2e, where we chose the negative root in order to obtain positive radial 
displacements (outward buckling). We also note that the maximum radial displacement 
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Pmax and the elastic energy are 




= —\J Scjly/ Ax 1.026 \/ Ax , and 

OTT 


e = 


^ J + 12^) di = Ai M />"' Ai . 


(35) 

(36) 


We finally discuss the existence of periodic solutions consisting of multiple identical folds 
of equal length. As before, each fold is symmetric (asymmetric in V’(5)) and its boundary 


conditions are 1) = = 0. Eq. (28) thus remains valid. The curvature at 


the boundary, however, is not zero anymore for periodic solutions (periodic solutions with 
Cl = 0 are unstable, as we will observe numerically in the next section). Denoting its value 
at s = 1 with Cl, the curvature boundary condition becomes 


y^cos - -^^sin f cot 


^ = 0 . 
0-2 


(37) 


To conserve the Hamiltonian, Eq (30) now yields the condition 


(38) 


Equations (37) and (38) leave the ratio Ci/a 2 of amplitude 02 and curvature Ci undetermined. 
We now proceed as follows: Since ci > 0 for radially outward buckling, and 02 < 0, we 


prescribe a negative ratio Ci/a 2 and search for roots of Eqns. (37, 38). Since we 


are only interested in the mode with lowest buckling force, we numerically determine roots 
of (37) and (38) in the vicinity of the zero-curvature solution, Eq. (32). The inset of Fig. 
shows that the minimal buckling force increases for Ci/a 2 < 0, and decreases for positive 
Ci/a 2 (dotted line). Since positive values of Ci/a 2 correspond to unphysical inward buckling, 
the inset suggests that the mode with lowest buckling force is the single fold with Ci/a 2 = 0, 
and periodic solutions are not selected at the buckling onset. This result is in contrast to 
the periodic undulations found for an axially compressed cylinder in absence of a substrate 


We can easily verify that the zero-curvature solutions with ci = 0 is a saddle point (or 
local extremum) of the buckling force by considering small perturbations around Ci = 0: 


Equations (37) and (38) geometrically describe two surfaces and S 2 in the space of 
independent parameters Ci/a 2 ). Both equations are simultaneously satisfied along 

the curve of intersection 7 of iSi and S 2 - Denoting the parameters of the single fold solution 
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as p* = (20/^3, (3/4)'/' TT, 0 ), the tangent T|p^ to 7 at p* is given by the line of intersection 
between the two tangent planes of and 32 - Basic geometry then yields 


T|p. = VSi\p, X V52|p, 


0,-(37r), 


~7r) 


(39) 


where V 5 i, 2 |p* are the gradients of 5 i ,2 evalnated at p* with respect to the parameters 
(/(■». s'"> , 01 / 02 ). Eq. ( |39| shows that the bnckling force does not change near p*, which 
agrees with having a saddle-point (or local extremnm) at p*. 


IV. NUMERICAL SOLUTIONS AND EXPERIMENTS 


To hnd the shape of the sheet for high deformations we solve the Hamilton eqnations ([^ 
nnmerically nsing a standard shooting method [ 2 Z]: for a hxed scaled force / and a trial 
valne for Pp = —fn at s = 0 the eqnations are integrated with a fonrth-order Rnnge Kntta 
method. For any trial fn, the valnes of V’, P, and p^ at s = 0 are obtained from the bonndary 
conditions (|^. The integration is stopped as soon as p = 0 is reached again. Every /„ which 
resnlts in 'ip{sc) = 0 corresponds to the prohle of a single fold (see Fig. 3). 

For low valnes of compression (Ax -C 1) the nonlinear solntions are symmetric and 


coincide with the ones of the linear regime as expected (see Fig. 3(a)). The associated 
compressive force / decreases with Ax. The negative sign of implies that solntions 
with more than one fold are nnstable: in eqnilibrinm the force / has to be a constant 
along the whole prohle for a hxed Ax. If the prohle consisted of N folds, each fold i wonld 
correspond to a solntion of the nonlinear system with (Ax)* = Ax/N {i G {!,..., A}). 
A small compressive hnctnation (Ax)* —?• (Ax)* -|- 6 wonld decrease the force /* needed 
to stabilize the fold i. Since the sheet is nnstretchable in the longitndinal direction, the 
adjacent fold wonld be less compressed than before. To keep it in place, however, a force 
which is higher than / and /* would be needed. Since this is not possible, the entire system 
becomes unstable. In particular, periodic solutions are unstable irrespective of the curvature 
at the boundary. Thus, unlike most other setups involving compressed sheets on deformable 
substrates HMZl [201EP, our system does not exhibit a transition from periodic wrinkles 
to localized folds; instead, it forms a single fold as soon as Ax 7 ^ 0. 

Above a critical axial displacement Axcrit ~ 1.88 the symmetry of the solutions is broken: 
the global energy minimum corresponds to an asymmetric fold. The associated force / 
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Aa: 
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(a) (b) 

FIG. 3: (a) Total elastic energy e (top) and horizontal force / (bottom) as a function of the 
compressive axial displacement Ax for the scaled nonlinear system (Eqns. 0/(©)- For small 
compressions the solutions are symmetric (blue solid branch) and coincide with the solutions of 
the linearized system (dashed black line (top) and black square (bottom)). Below a compressive 
force of approx. 7.8, asymmetric solutions appear (red branch) which are energetically favoured 
over the symmetric ones (blue dotted branch). The two branches are depicted up to the point 
of self-contact (assuming a vanishing thickness h = 0 of the sheet). Insets: Profiles for different 
values of Ax. (b) Normal force at the circular boundaries, fn, as a function of the horizontal force 
/ for the scaled nonlinear system (Eqns. 0/(i)- Its value is the same on both boundaries when 
the fold is symmetric. For the antisymmetric profile two different values are found when fixing /. 
The higher one corresponds to the boundary to which the fold has tipped over. 


decreases quicker than before and changes sign at Ax ~ 2.22. At this displacement value 
the energy displays a maximum and the corresponding external force is zero. For even 
higher Ax we obtain a recumbent fold and one has to pull instead of to compress to hold the 
asymmetric sheet in equilibrium. When the sheet is confined further it starts to touch itself. 
Assuming a vanishing thickness h, one finds Ax ~ 2.62 for the asymmetric and Ax ~ 2.54 
for the symmetric branch. 

We validate our theoretical results (f) with the help of finite element simulations (see 
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(a) (b) 

FIG. 4: (a) Maximum radial displacement Pmax as a function of Ax for a neoprene sheet jackating 
a cylinder of radius Rq (dots). To compare the results to the nonlinear solution p{s = Sc/2) of the 
scaled system (solid blue line), all measurements were divided by \/0.2cm Rq. The black dashed 
line shows the solution in the linear regime (Eq. ®). (b) Comparison between experimental 
(points) and theoretical nonlinear profiles (dotted lines) for Ro = 3.75cm and different values of 
AX. All profiles are plotted in unsealed units. The thickness h = 0.2cm of the sheet is indicated 
with double arrows. 


appendix), and {ii) in a series of experiments with neoprene sheets wrapped around cylinders 
of different radii. In these experiments, the sheet is prepared in such a way that it is free of 
stretching. The experimental prohles are then photographed and extracted with Imaged |38] . 
While relatively simple, we note that this technique only allows the analysis of symmetric 
prohles before the tipping point. [H] Fig. shows the experimental results together with 
the theoretical predictions. Since all material parameters scale out of the theory, no htting 
parameters are needed. In Fig. |4(a)| the maximum value of the radial displacement, pmax, is 
plotted as a function of the axial compression Ax. One observes that the experimental curves 
converge towards the theoretical solution for increasing radius Rq. This can be explained 
by taking a look at our theoretical assumptions again: to simplify the problem we have 
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assumed that -Ro/h is large (see Sec. The thickness of the neoprene sheet is constant and 
approximately h = 0.2cm in all experiments, implying that our theoretical approximations 
become more accurate for larger radii Rq. The effect of the cylinder radius is also noticeable 
when looking at the tipping point: while theoretically predicted to occur at Ax ~ 1.88, 
experiments yield smaller values that decrease with Rq. 


Fig. 4(b) shows three different experimental prohles together with the corresponding 
nonlinear theoretical solutions for a cylinder radius of Rq = 3.75cm. Note that experimental 
and theoretical results are shifted radially, since the digitized prohles describe the outer 
side of the sheet {h = 0.2cm away from the solid cylinder) whereas the nonlinear solutions 
correspond to the theoretical centerline of the sheet. Despite the simplihed nature of our 
model, theory and experiment coincide remarkably well even for high curvatures at the tip 
of the fold. 


V. CONCLUSION 

Euler buckling is a well-studied phenomenon with numerous applications and occurs 
in various situations in nature and engineering |3l]. In its simplest form it deals with the 
situation of translatory invariance in one direction, such as the buckling on planar substrates 
[T^ . However, buckling often occurs in curved geometries inducing non-isometric 

deformations of the sheet. In this article we have exemplihed the effect of curvature on 
wrinkling by considering an Euler-type buckling in a cylindrical geometry. In our setup, an 
elastic membrane is wrapped around a cylindrical substrate of same radius and compressed 
axially. Instead of a wrinkle-to-fold transition typical for elastic substrates, a single fold 
appears immediately above the buckling threshold. For even larger compression, the fold 
breaks symmetry, leading eventually to negative buckling forces. The theoretical predictions 
of our simple model were validated with experiments on neoprene sheets and hnite element 
simulations. Without a htting parameter, experiments, simulations and theory coincide 
rather remarkably. 

While we do not observe wrinkling solutions, we speculate that such periodic solutions 
would form only if the radius of the unstretched elastic sheet was larger than Rq\ For 
small amplitudes, the sheet would not be in contact with the cylindrical substrate and 
thus form the periodic undulations of a free, axially compressed cylindrical sheet. For larger 
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compression, contact with the rigid substrate eventually occurs, which then most likely leads 
to the single fold solutions discussed in this article. 

In the experiments one can force the sheet over the point of self-contact. Due to the 
friction between the sheet and the substrate a stable equilibrium can be found whose Ax 
depends on the value of the friction coefficient. Our system can thus be interpreted as a 
simple mechanical switch with two stable states: the cylindrical geometry and an asymmetric 
conhguration with Ax > 2.62. To switch from the cylindrical state to the other, an external 
compressive force above the critical buckling force / = 11.55 has to be applied. To return to 
the unstretched sheet requires stretching the fold (decreasing Ax) into a regime of positive 
/ > 0, after which it will unfold on its own. Consisting only of a single movable part 
and featuring simple assembly steps, such a switch could potentially be used in micro- and 
nanoscale applications. 
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Appendix A 

1. Derivation of the elastic energy 

We derive the elastic energy of the cylindrical shell by hrst linearizing the elastic strains 
associated with a deformation of the cylinder. Comparing the deformed and undeformed 
inhnitesimal length of two neighbouring material points, one hnds for the stretching strains 
at linear order in R: 


= 0 , 

R — Rq 


^ipip — 


Rn 


(Al) 

(A2) 


where (p measures the angle in the circumferential direction on the cylinder. Note that the 
hrst equation follows from the inextensibility in the F-direction (see next section), and the 
off-diagonal es^p = 0 for axisymmetric twistless shells. Similarly one obtains for the bending 
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strains 


Kss = , 




= 


COS V’ 
R 


1 

Rq 


(A3) 

(A4) 


where and are the curvatures in the circumferential direction of the deformed and 
undeformed conhguration, respectively. Due to axisymmetry, the off-diagonal strain again 
vanishes, Ks,p = 0. Using Eq. (A2), we write 1 


l + Ciz 


Rq 


. Since the bending 

strains will only dominate for deformations close to isometric ones, where the stretching 
strains are ~ 0, we may further approximate ~ ^ [3l], leading to the hnal 

expressions for the linearized strains 



R — Rq 

(A5) 

Ro ’ 

Kss = 


(A6) 


cos'll) 1 

Ro Ro 

(A7) 

RS(p 

ess = esv? = 0 • 

(A8) 


Using an isotropic Hookean material with Young’s modulus Y and Poisson ratio z/, the 
simplest form of energy densities with decoupled bending and stretching contributions reads 


2('x _ i^ss + + 2z/esse^<^) , 

Yh^ 

24(1 - z/2) 2 uKssK^^) . 


(A9) 

(AlO) 


The total energy is then obtained by integration over the undeformed reference surface. 


E 




4>(e5S, ^(A'ss, K^^)Rq dtp dS 


from which Eq. Q is obtained. 


(All) 


2. Comparison with finite element simulations 

The mechanical model leading to Eq. ([^ assumes linearity of strains, inextensibility along 
the axial direction, and radii of curvatures that are large with respect to the shell thickness h. 
To qualitatively test the validity of these assumptions, we compare the deformed membrane 
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(a) (b) 

FIG. 5: (a) Abaqus simulation of a recumbent fold for Rq = 1.0, h = 0.01 and AX = 0.232. 

Color-coded is the relative radial displacement {R — Ro)/Ro. Only a segment of the cylindrical 
membrane is shown for better clarity, (b) Comparison between Abaqus profiles (points) and 
theoretical nonlinear profiles (dotted lines) for two values of Ax. All profiles are plotted in scaled 
units. The thickness of the sheet is indicated with double arrows. 

profiles of Eq. ([^ with numerical simulations obtained from the commercially available 
Abaqus hnite element software package [SH]. In Abaqus, the membrane is modeled as a 
three-dimensional body using the reduced-integration quadrilateral elements CAX4R, which 
take the axisymmetry of the problem into account to reduce computational cost. We found 
that a spatial resolution of 2 elements in the thickness direction is sufficient to capture the 
deformation accurately, with higher resolutions leading to no observable change in the prohle 
geometry. The rigid cylindrical constraint is modeled using Hertzian contact dynamics, 
penalizing hnite element nodes that impenetrate the cylinder |1D]. Since the fold length 
Sc changes during axial compression but is not known a priori, we simulated a cylindrical 
membrane of length L 3> Sc- The geometries considered are h/R^ between 5 ■ 10“^ and 
10“^. At the membrane boundary, i?(0) = R{L) = Rq and ^’(0) = t/’(T) = 0 is imposed. 
To compress the membrane, we prescribe the axial displacement Ux by enforcing 14(0) = 0 
and Ux{L) = —AX. Starting from AX = 0, we slowly increase AX to obtain fold prohles 
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at various stages of compression. 


Fig. 5(a) shows an axisymmetric fold obtained from a typical simulation (for better clarity, 
only a part of the membrane is shown), agreeing qualitatively with the recumbent fold 
found in experiments and using the simplihed model (|^. To investigate the validity of axial 
inextensibility, we compare prohles in the early and intermediate stages of folding, where 
errors introduced by inextensibility are expected to be maximal due to the large compressive 
forces present in the fold. Fig. |5(b) shows that our simplihed model compares well to the 
prohles of Abaqus even for large compressive forces, demonstrating, in particular, that 
inextensibility is a valid model assumption. 
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